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Abstract 

Thermoelastic phenomena associated with orientational phase transition in 
solid Ceo are considered. Coupling of the order parameter with elastic strain 
is analyzed in the Landau theory of phase transitions. Wide range of possible 
coexistence of the FCC phase with SC one in Cgo has been found from 260 
K down to critical temperature 160 K of the instability of high-symmetry 
phase. The specific heat, thermal expansion coefficient and isothermal com- 
pressibility for low-symmetry phase are calculated as function of temperature 
and pressure in an agreement with reported data. Applicability of the Landau 

theory to ordering transition in fullerenes is discussed. 
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Among other crystalline fullerenes solid Ceo exhibits orientational ordering transition [|TJ . 
There are extensive experimental data available about crystallography of this transition, 
associated anomalies of bulk modulus and thermal expansion coefficient as well as effect of 
applied external pressure, however, no adequate theory of such thermoelastic phenomena has 
been developed so far. Phenomenological approach based on the mean-field Landau theory 
of ordering transition is used in the present paper to analyze these effects. Due to high 
symmetry of the buckyball Ceo molecules microscopical theory of ordering transition involves 
symmetry adapted rotator functions to describe orientational interaction and fluctuations 
@. However, spontaneous lattice distortion originates from coupling of the orientational 
order parameter with homogeneous strain and, hence, mean-field Landau theory ||||] can 
be employed. 

Thermodynamic potential around the symmetry breaking phase transition can be ex- 
panded in this theory || in powers of 'order parameter' rj. The theory initially was de- 
veloped for second-order transition where order parameter is continuous at the transition 
point, however weakly-discontinuous first-order transitions like the orientational ordering 
transition in solid Ceo can be considered as well || . If the symmetry groups of both parent 
and product phases are known a priory then scalar order parameter can be used and the 
Ginzburg-Landau expansion of the free energy difference between the parent and product 
phases has the form |5]|| 

A£ = |(T - T C ) V 2 + |r? 3 + ^ 4 - Ei]. (1) 

where T c is a critical temperature and E is external conjugated field. Only second-degree 
coefficient is supposed to depend on temperature and equilibrium value of r] is determined 
by the minimization of AQ. The stability requires highest order coefficient C to be positive 
and the third degree term B^O implies first order of the transition. 

In the E = case the Gibbs free energy (||) has two possible minima. One with r] = 
corresponds to a high-temperature undistorted phase stable for T > T c . For 



low-symmetry phase exists with 

B I /Tn-T\5 
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The phase energies become equal at the temperature of first-order transition T* = T c + |^ 
where 77 jumps to —\%- The temperature range between T c and To where both high- and 
low symmetry phase are stable and can, thus, co-exist is an absolute thermal hysteresis, 
real hysteresis being determined through the nucleation and growth kinetics by actual cool- 
ing/heating rate. 

The temperature dependence of specific heat in the low-symmetry phase 

AC P (T) = -T d ~ M ' 



U^( T ^Y\ (3) 



g T 2 

2 C ( X ' V T - T 

implies divergence oc (T — T)~^ as T — > To. This can be seen as some sort of 'precursor' 
phenomena below the transition temperature in a contrast with the critical region near the 
continuous transition, where specific heat anomaly due to critical fluctuations around T c 

<" 2 > - • (4) 

appears both in high- and in low-symmetry state. 

Full review of the Landau theory of ordering transitions in crystalline fullerenes can be 
found elsewhere ||. The transition in Ceo takes place around 260 K and reduces crys- 
tallographic symmetry from FCC (or Fm3m) lattice into SC one with space group Pa3 
|7],|8|]. The phenomenological order parameter rj defined by the symmetry change belongs to 
representation of Fm3m space group corresponding to X point of the Brillouin zone 0]. 

To analyze thermoelastic phenomena the coupling of 77 with homogeneous strain tensor e 
should be considered and the Ginzburg-Landau expansion should include the terms associ- 
ated with elastic energy. It was shown that symmetry allows volumetric strain eo = Tr(e) 
to be coupled with rj in the lowest order by term 
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AQ int = DotoV 2 ■ (5) 

Any shear strain has a higher degree of coupling with the order parameter and can be ne- 
glected in the present consideration. Furthermore, there might be no external homogeneous 
field conjugated to the order parameter associated with wave vector q ^ 0. Coupling be- 
tween the order parameter and volume change similar to Eq. fl5|) has been recently considered 
in phenomenological theory of weakly-discontinuous ferroelastic phase transition in cu- 
bic crystals. It was shown that the minimization of free energy with respect to eo leads to 
renormed Ginzburg-Landau expansion of AQ in powers of rj ||: 

a' 

Ag(T,P, V ) = Ag + -(T-T>)r ] 2 

+jV 3 + ^V 4 ~E V , (6) 

with coefficients 

a' = a + 2a v D (7a) 
9 n 

K = T C + ^P (7b) 
C' = C-^. (7c) 



Here AQ (T, P) is thermal expansion energy of a high-symmetry phase with r] = [|1C 



A£ (e ) = -a v K (T - T c ) e + ^ + Pe , (8) 

with a v and Kq being the volume thermal expansion coefficient and bulk modulus of the 
high-symmetry phase, respectively. The applied hydrostatic pressure P shifts the transition 
temperature: 



(9) 



dP K a> 

The expressions for transition-related anomalies in isothermal compressibility and ther- 
mal expansion coefficient can be obtained |9| by substituting equilibrium 77 (T, P) for low- 
symmetry phase into Ginzburg-Landau expansion and then differentiating with respect to 



temperature and/or pressure. The difference in thermal expansion coefficient between low- 
and high-symmetry phase has the following dependence on the temperature 

^v(T) = (AS - A£o) 

This expression along with analogous one for compressibility are similar to Eq.([3]) for 
ACp and imply 'precursor'-like divergence near To. The discontinuities at the temperature 
TL of the first-order transition 



AC -2^ Aa - A3 - ^ (11) 

AC P -2T—, &oc v -—, Wt-c^i (11) 



are related by the Keesom-Ehrenfest relationships || 

dT' c A(3 T TlAa x 



dP Aa v AC 



(12) 



Several different values for ^ between 105 and 160 K/GPa were obtained in different 
studies of Ceo with most of the data located around 110 K/GPa |iTHi"5||. Bulk modulus of 
the FCC phase was reported to be 16.4 GPa [14]] that is intermediate between other values 



Pqjlq| . For low-symmetry phase the bulk modulus was found to be 11.6 GPa p?7| . It gives 
the discontinuity of isothermal compressibility A(3t ~ 2.5xl0 -2 GPa -1 and implies Aa v ~ 
22.5xl0 -5 K _1 with ACp ~ 0.23 kJ/(mol K) through the Keesom-Ehrenfest relationship 
Eq . (|l~2"D . This Aa v is in qualitative accord with results of the diffraction study |2(| where 
a v for FCC phase was found to be 6.1 xlO -5 K -1 and definite precursor phenomena were 
observed in low-temperature phase with average a v ~ 21±3xl0~ 5 K -1 . Our value of ACp 



appears to be twice as much as experimental value [JT8 . 

The dimensionless volume effect Aeo was found around -0.01 [|l~9| , f2"0| and it gives the 
absolute hysteresis 

^ ^ B 2 9 Ae dT c 

T n — T = = w 99 K 

c Aa'C 4A/3 T dP yJ 



Taking the highest reported transition temperature 261.5 K [2~I| as Tq we get T' c w 162.5 



K and T^(P = 0) ~ 251 K. This is in very good agreement with DSC experiments [22 



where two peaks around 250 K and 150 K were found in Cp of Cqq upon cooling and only 
sharp peak at 260 K has been found during heating. Peak at 150 K can be attributed to the 
critical order parameter fluctuations in supercooled high-temperature phase near T c whereas 
peak around 250 K corresponds to the onset of the transition accompanied by increase in 
Cp of low-symmetry phase near T . This peak increases in superheated SC phase during 
heating cycle. Sharp peak in the sound velocity and attenuation was found near 160 K 
in Ceo single-crystal and the anomaly in the X-ray diffraction peak intensity was also 
reported around 155 K [Q. 

The temperature dependence of a v in low-temperature phase plotted according to 
Eq.(|10"D is shown in Fig.p], where qualitative agreement with the experimental data |25 



can be seen. Similar expression appears for specific heat and isothermal compressibility, and 
the temperature dependence of bulk modulus in SC phase shown in Fig.|2] can be obtained 
accordingly. Pressure acts through the change of T' c by Eq. ( [7b]) and the pressure dependence 
of bulk modulus in the low-symmetry SC phase 

*(t,p)=4i+||(i+ 



r--%!£gpy „ , (13) 



\T -T C K B 2 

is plotted in Fig.|3] for room temperature. Eq. (|13"D is essentially different from the Murnaghan 
equation K(P) = K(0) + K'P used in some experimental studies to describe measured data. 
The assumption of linear dependence of bulk modulus on pressure has led to conclusion about 
negative jump in compressibility upon transition from FCC to low-symmetry SC phase at 



P = |26] even though in the same study positive compressibility jump was found at 293 K 
for observed pressure-induced transition around 0.2 GPa. As the transition discontinuities 
(|TT|) do not depend on applied hydrostatic pressure and volume effect has been indeed found 
independent on P in that experiment, one can suggest that increase in K is a result of linear 
extrapolation of obtained data outside the region of the stability of SC phase at P = 0. 
The Murnaghan equation was found inadequate in recent pressure experiment [^7] as well, 
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where bulk modulus of SC phase at 170 K was found to be 14.2 GPa, that is in qualitative 
agreement with our K(P) curve in Fig.(|3|). 

Thus, thermoelastic singularities associated with the orientational ordering transition in 
crystalline fullerene Ceo can be reproduced in the Landau theory through temperature and 
pressure dependence of the Gibbs free energy in SC phase and coupling of elastic strain with 
order parameter. The Keesom-Ehrenfest relations (|l~2[) for second derivatives of the Gibbs 
free energy appear to be satisfied, hence, the energy is continuous around transition and 
can, thus, be naturally written in the Ginzburg-Landau form due to symmetry breaking. 

Usual belief is that the singularities around phase transition are due to effect of critical 
fluctuations and for many transitions of the 'order-disorder' kind the fluctuations are known 
to be strong. Thus, agreement obtained in the present approach within the Landau theory 
which ignore the fluctuations implies question why they do not make significant contribution 
in this case. Possible answer is that the transition temperature Tl is located far from critical 
temperature T c , where high-temperature phase becomes unstable and considerable fluctua- 
tions arise. One can suppose that the coupling of the order parameter with strain suppress 
inhomogeneous fluctuations because elastic forces have long-range behavior and the strain 
energy associated with order parameter fluctuations is large. In order to make quantita- 
tive estimation the experimental information concerning correlation length for orientational 
fluctuations is needed. 

We can evaluate, however, the homogeneous order parameter fluctuations in the present 
model. In order the weakly-first-order transition to be properly identified and clearly sep- 
arated from the background of fluctuations in high-temperature phase, the discontinuity of 
the order parameter should be greater than its thermal fluctuations It gives us criterion 
for the third-order coefficient: 



B > 




(14) 



We have 



B*_ 

C' 3 



81 (Ae ) 2 
2 Af3 T 



69.6 kJ/mol , 
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that is much larger than the energy of thermal fluctuations ksT^ 2.08 kJ/mol and, hence, 
the homogeneous order parameter fluctuations are very weak near T' m . 

To conclude, we have analyzed the thermoelastic phenomena around the orientational 
ordering transition in crystalline fullerene Ceo m the frame of the mean-field Landau theory. 
The transition discontinuities of thermal expansion coefficient and isothermal compressibility 
originate from the temperature and pressure dependence of the Gibbs free energy in low- 
symmetry phase through coupling of elastic homogeneous strain with order parameter. The 
temperature range of the coexistence of FCC and SC phase appears to be very large - from 
260 K down to critical temperature T c pa 150 K and it can explain experimentally observed 
anomalies around this point. 
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FIGURES 
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FIG. 1. Temperature dependence of the volume thermal expansion coefficient in low-symmetry 
phase. 
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FIG. 3. Pressure dependence of bulk modulus in SC phase for room temperature 25°C. 
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